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COOLING AND RECOMBINATION PROCESSES IN COMETARY PLASMA 

M. K. Wallis and R. S. B. Ong 

1. INTRODUCTION 

It has long been recognized that collisional 

cooling and recombination processes are likely to be 

4 

important in the inner cometary coma, in a 10 km radius 
region for the larger comets (Biermann & Trefftz 1964) 
Cherednichenko (1970) laid stress on dissociative 
recombination processes, as possibly playing a role in 
the production of observed ions and radicals. Oppenheimer 
in his spirited contribution to this conference, 
emphasized that a variety of ion-molecule interactions 
occur relatively rapidly and probably take part in the 
production of known cometary radicals. 

In this paper, we focus our attention on the ion- 
electron plasma in comets and examine in the first place 
the cooling processes which result from its interactions 
with the neutral coma. For the plasma is generally 
very energetic (1-lOOeV) and must be cooled if it is to 
reach moderate densities and promote efficient particle- 
particle interactions. For example, solar wind electrons 
have 10-15eV energy, they experience some adiabatic 
heating (factor 2 or 3) in passing through the coma, they 
may gain around lOeV in passing through a collision- 
free, resistive shock and perhaps suffer additional 
heating via plasma turbulence effects. Photo-ionization 
processes may release other energetic electrons — He 
584^ photons could give electrons with about lOeV 
(Biermann & Trefftz 1964), although most have less than 5eV 


856 


I 


5 6 

New cometary ions produced at 10 -10 km in the far 
coma probably gain most of the streaming energy of the 
solar wind, through being accelerated in the E and B 
fields up to perhaps several thousand eV (Wallis 1973a). 
How quickly these ions are lost from the incoming solar 
wind plasma largely determines the ion pressure. 

Cooling processes have general relevance for plasma 
behaviour in comets , in describing the overall plasma 
flow through the coma and in cometary plasma formation. 
Specific problems that have received attention and 
require a careful description of the cooling rate are 
that the visible ion structures cannot consist of hot 
and therefore low density plasma; that cool molecular- 
ion plasma is rapidly destroyed by dissociative 
recombination;’ and that energetic photo-electrons 
would exert a high pressure in the inner coma and prevent 
penetration by the solar wind. We develop a continuous 
description of the cooling effects in order to look at 
such problems. 

In this preliminary examination , we shall consider 

a cometary coma composed predominantly of H 2 O and its 

decomposition products (Wallis 1973b). For specific 

estimates, we use a comet of the size of comet Bennett 

9 9 

1970 II, with a production rate Q = 10 H 2 O molecules 
- 1-1 

ster s at 0.7a.u. heliocentric distance. The coma 
density depends a little on assumptions about the 
expansion velocity V ; this factor is relatively 
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unimportant, but for concreteness and consistency, we 
suppose V increases with distance due to photo- 
dissociative heating (Wallis 1974), so that the density 
at radius R is 


N 


= _SL 

VR^ 


where 
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at R = 10^-10*^)011, 


COOLING PROCESSES 


Descriptions of electron cooling are given in planet 
ionosphere studies (Henry & McElroy 1968, Sawada et al. 
1972, Olivero et al . 1972), energy loss rates in 0, CO, H^O 
etc. being computed on a continuous slowing-down 
approximation. Data for e-OH collisions are incomplete* 
and we suppose it comparable to CO above 7eV , while 
similar to H 2 O with rotational transitions dominant at 
lower energies (Shimizu 1974). Ionization data for 
H 2 O, OH and 0 have been summarized by Wallis (1973b). 

Solar wind protons and energetic cometary ions are 
lost from the plasma primarily in charge exchange 
processes with neutral gas , having cross-sections 
a = l-3xl0"'^^cm^ at 10^-10*^eV. 

Electrons are cooled in a variety of processes at 
rates varying with energy as shown in Fig.l. The 
functions shown are a continuous approximation to the 
discrete energy losses actually occurring which is useful 
in calculations (Olivero et al.l972). The approximation 

exaggerates the width of the 'holes' in the CO cooling 
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Fig. 1 . Cooling rates for inelastic collisions of 
electrons with energy e in molecular gases , on the 
continuous slowing-down approjcimation . The theoretical 
cooling rate via rotational excitations of OH is similar 
to that of H 2 O (section to the left), but both are 
uncertain to a factor 3. The structured, low energy 
part of the CO curve is due to vibrational excitations. 




function, where this becomes small or even zero. But 
after making allowance for the redistribution of 
energy between the electrons , which has to be done 
anyway, the inaccuracy due to the continuous 
approximation becomes small. For electron collisions 
with H^O, rotational excitations dominate below 5eV 
and as each energy jump is small, the continuous 
approximation is good even for single electrons. 

The cooling rate , calculated on the rigid rotor 
approximation for electrons of energy e exceeding 
A.e = 0.0 2 5eV has the form (personal communication from 
M. Shimizu) 

de/dt = No .y Ae = -as , a == 5x10 ^eV^^^cm^s 
rot e ’ 

Above 5eV, electronic excitations and ionizations become 
important (Fig.l) and the cooling rate increases 
steeply in 6-20eV as 

de/dt = -a"Ce-3eV)^N , a' = 2xl0~^ eV'^cm^s"^. (3) 


The time for cooling to the minimum energy Ae depends 
little on the initial energy if above lOeV: 
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For electrons in CO, the cooling function is 
significantly structured (Fig.l), particularly because 
of the sharply-peaked vibrational excitations below 
5eV. It is more meaningful to calculate the average 
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cooling rate over a Boltzmann distribution, which turns 
out to be approximately linear above ^eV : 

dF/dt == -bFN , b = 1-1.5x10 ^cm^s (5) 

Expression (5) is only applicable if thermalization 
processes are rapid enough. The thermalization rate 

due to Coulomb collisions is 

^ -1 ^-3/2 Q^-in"5 ,,3/2 3 -1 

t = cne , c = 8x10 eV cm s , (.6) 

c ’ 

which has to exceed the cooling rate of (5): 

t ^ > bN or n/N > (7) 

c 

In practice, condition (7) is not, fulfilled for 5eV 
electrons at densities found in the inner coma (Table 1). 
Plasma instabilities will therefore play a role in 
thermalization. For a highly anisotropic velocity 
distribution , the thermalization rate is a fraction of 
the plasma frequency (Davidson 1972) 

-1 1 8/2 -1 
t =0.1u) =dn^, d = lOcm'^ '^s . (8) 

~ pe ’ 

. . -1 

For the densities of Table 1, t^ exceeds bN by more 
than a factor 10. A more detailed treatment would 
modify (8) to include the damping effect of electron- 
molecule collisions, but this is estimated to be 

3 

significant only inside 10 km radius. The conclusion 
is that the anisotropy in electron energies will 
generate plasma turbulence, which produces some 
thermalization and limits the growth of anisotropy: 
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Table 1 

Major ionization and recombination paths. 

E2O 


OH 








r ' d 

5xl0"^° 

5xl0~^° 

10 

3 -1 

3 -1 

6 -1 

cm s ^ 

cm s 

cm s 



r 

H 01 

i+H H. 

,0+H 2 H 2 O 




0 . 5-lxl0“^e"^ 


^ 5.5xl0"'^e'^ 


— 7 — 1 

5 2 . 2 ±0.4x10 e ^ 


4xl0"'^e~ 


The data is culled from Banks & Kockarts (1973), Cherednichenko (1970) and 
Leu et al. (1973). For electron energies exceeding 0. leV, the dissociative 
recombination coefficients a may decrease as sharply as e~2. At the relevant 
cometary densities, the 3-body collisions producing H 2 O . nHgO are negligible. 
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expression (5) should remain an adequate approximation 
to the cooling rate. 

The electron cooling functions (.2) and (3) do not 
tend to give large anisotropies. Plasma instabilities 
are less likely to be important in thermalization and 
the unaveraged cooling rates should be appropriate. 

3. IONOSPHERE OF THE H 2 O COMET. 

Suppose conditions in the inner coma are quasi- 
stationary with photo-ionizations being balanced by 
recombinations , changes due to the outward expansion 
being comparatively slow. Photo-ionization releases a 
spectrum of electrons, mainly below 5eV, which cool 
through vibrational excitations of the H^O as (2) and 
subsequently recombine dissociatively . The ions may 
undergo ion-molecule interactions before recombining 
as shown in Table 1, but this makes little difference to 
estimates of plasma density. 

The H 2 O photo-ionization cross-section (Metzger & 
Cook 1964, Wantanabe & Jursa 1964) varies little over 
12.6-18.3eV, so the electron energy spectrum in 0-5. 7eV 
is close to the (shifted) solar spectrum: 

P(e) Nt’^ j (1 + e/e.,.)”^”^, e... = 12.6eV. (9) 

Here we have used a power law approximation to the solar 
spectrum (e.g. Shul'man 1972) with index j = 8.7. With 

the continuous description of the cooling (2), the 
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electron distribution function for energies above 
thermal where recombination is negligible is 

f(e) = ^ PU')dc' = -4 e^Cl + — ) ^ . (10) 

aN ax 

e 

o 

The solar 584A photons contribute a further 10% by 

number to PCe) probably mainly in dissociative ionizations 

releasing electrons of 2-3eV energy, so are relatively 

o 

unimportant. The 30 4A and other far UV photons 

contribute fewer and more energetic electrons , which 

are very rapidly cooled according to (3) so are also 

unimportant. It is thus adequate to use expression (10) 

for the hot electrons , whose total number density and 

pressure do not exceed 

, 3/2 tt/2 

f ^ it f 2k-4 3 

f(e) de cos^ ^0d0/ax = 40/cm 

Jq j - 1 J 0 

( 11 ) 

% r~ Q 

ef(e) de _< f(e) de = 35 eV/cm . 

0 2j-3 -'O 

The numerical estimates apply for the ionization time 
6 

X = 10 s (at O.Va.u). They are uncertain by a factor 
2 or more because of uncertainty in the constant a, 
taken from C2 ) . 

Let us suppose for simplicity that recombination 
processes occur relatively slowly, so that recombination 
occurs only subsequent to cooling to thermal energies. 

This is valid if the density of thermal electrons is far 
greater than the density of energetic ones by (11). 
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Neglecting transport effects, the steady state balance 
of ionization and recombination rates is then 

N^t = an^ 

where a is the relevant dissociative recombination 
coefficient of Table 1. Thus the density of thermal 
electrons is 

n = (N/aT)2 = 0.5-1. 5 cm^^^. 

Specific values of nCfor a = 2.5x10 As ^cm eV^s ) are 
given in Table 2, and clearly exceed the density (11) 
of the energetic electrons inside lO^km radius. The 
lifetime of an ion before recombination is (an) = i(n/N), 
which is far shorter than an e-folding time for changes 
due to the flow, R/2.3V, so ion transport effects are 
indeed negligible. It can also be confirmed that electron 
thermal conductivity is adequately limited by e-H^O 
collisions . 

In the inner coma, this ionospheric plasma is 
closely coupled to the neutral gas and streams radially 
outwards with it. Outside some radius of the order of 

Qa/V 3xl0^km, (13) 

ion-molecule collisions become infrequent and the plasma 
can behave as a separate fluid with a smaller mean free 
path fixed by gyro-radius or collective plasma effects. 

A tangential discontinuity might exist between the 
ionospheric plasma and the plasma of solar wind origin 
(Wallis 1973b, Schmidt 1974). The ionospheric ions and 
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electrons would exert a pressure of the order of their 
stagnation pressure 

t f(e) de + nCAe + kT^ + ) . (14) 

The first term representing the suprathermal electrons 
is given by (11) and, on equating to the neutral gas 
temperature T, the values of P(stag) given in Table 2 are 
found. 

This quantity P(stag) is to be compared with the 

solar wind stagnation pressure, which is of order 
4 - 3 

10 eV cm at 0.7a.u. This would place the ionosphere 

3 

discontinuity and stagnation flow region within 10 km 
sunward of the comet, impossibly far inside the 
decoupling radius (13). If the dissociative heating and 
and enhanced expansion velocity of the present model coma 
are discounted, the values of P(stag) would be lower and 
the concept of a plasma contact discontinuity still more 
dubious . 

^ 4. PLASMA FLOW THROUGH THE COMET COMA . 

As long as the ionospheric plasma pressure is low 
(section 3), the solar wind plasma can flow on into the 
inner coma. We consider that it picks up new cometary 
ions and loses those neutralized in charge exchange 
processes in interactions with the neutral coma 
(Wallis 1973a). There is no sudden change at the scale 
coupling radius Qa/V (13), and we can expect the flow to 
penetrate far inside this position. We are interested 
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Table 2 


Parameters in an HgO coma. 


R 

(km) 

V 

(knu/ s ) 

T 

C 

kT+mV^ 

(eV) 

N 

( cm“^ ) 

n 

( cm~^ ) 

Pstag 
(eVy cm^ ) 

10^ 

0 .55 

350 

0.08 

2x10^° 

1.1x10®, 

1.2x10^ 

iq3 

0.8 

750 

0.18 

1.3x10® 

8.6x10® 

1. 8x10® 

10^ 

1.5 

2800 

0.67 

7x10® 

6. 3x10® 

4.7x10® 

o 

cn 

2 . 6 

C650 ) 

1.27 

4x10® 

50 

100 


Source strength Q = 1 q 29 molecules ster'^s* . The neutral gas expansion 
velocity V, density N and temperature T are taken from the heated coma model 
of Wallis (1974), but T at lO^m is uncertain. The plasma density n and 
stagnation pressure are calculated for solar radiation conditions at 0.7 a.u. 
heliocentric distance. 



here in how rapidly the plasma can cool and condense. 

There is no flow solution yet available for this 
strongly interacting and strongly cooled plasma flow, 
so we shall in making definite estimates assume that 
the plasma velocity in the incoming flow sunwards of 
the comet is 

u = R/t^, = lO^s, (15) 

being an empirical flow time scale which fits in 
with the outer flow solutions in the l-5xl0^km region 
CWallis 1973a, b). The assumption (15) is not 

critically important: if in^ error, the distance scale 

that we derive is simply distorted. 

The cometary ions are rather energetic: they take 

up most of the streaming energy of the decelerating 
flow and have gyration velocities of the order of the 
flow velocity v.,. = t : at the place where they become 

ions . They are lost primarily through charge exchange , 
so their distribution function gCv..,) satisfies 


^ = 
dt 


-u 


vr2 


V... 


g- 


(16) 


This equation neglects increases in ion energies due to 
continuing adiabatic compression. The solution to (16) 
using (15) and taking V as constant is 

g(v.,. ) 'V exp i2.^^(R,,. ^ - R ^), (17) 


where 


= ii-(v.,.) 
1 


= { 


o(v^) V.,. 
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The number is reduced by a factor e by the 
position R = and a further factor e by 

4 

R = 0.7 . Numerically xs 3x10 km for inxtial 

4 

ions of OOOkmys velocity and = 2x10 km for the 50km/s 

4 

ions formed at 5x10 km. The solar wind protons have 
thermal speeds of the order of 50-100km/s , so the 

4 

corresponding disappearance scale is = 2-2.5x10 km. 

The electrons in the inflowing plasma also cool 
rapidly due to various ionization and excitation processes. 
We suppose the cooling rates of (2) and (5) are 

4 5 

representative in the mainly 0 and OH coma in 10 -10 km. 

5 

The electron energy is perhaps £Q=50eV at Rq = 10 km 
and decreases as 

I/Eq = exp - ^ > 6eV, (18) 

according to (5) with C15), the scale radius being 

^e ■ ^"^f 2xlo\m. (19) 

4 

The mean energy reaches 6eV at R^ = 1.5x10 km and would 

4 

be leV at 10 km by formula (18), but even faster 
cooling according to expression (2) is appropriate: 

(6eV)^^^ - (R~^ - R, "^ ) . (20) 


4 V f 1 

The electrons become fully cooled, it follows, at the 
position 

{R^"^ + (6eV)^^^4V/3QaT^}"^ l.OxlO^km. (21) 

So in the absence of heating mechanisms , such as plasma 
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turbulence transferring energy from the ions, the 
electrons cool explosively fast between positions — 



If plasma is to flow from the coma out laterally 
into tail rays, it is clear that the same scale radii 
are important. For example, suppose that flow occurs 
at constant radius and speed Cthe pressure gradient 
balancing the effective friction). 

We replace d/dt in (16) by Rt d/ds and obtain 
g -v exp - Cs-Sq )/R^ . 

With flow distance s-s^ = nR/2 , we see that most of 
the ions would be lost if R < Similarly, the 

electrons would be strongly cooled if the lateral flow 
takes place at Coincidentally, these ion and 

electron scales are very similar in magnitude. 

5. DISCUSSION 

Solar plasma plus accumulated cometary ions and 

electrons is affected very strongly as it flows into 

4 4 

the coma from 2x10 to 10 km (this value for the comet 

29 -1 -1 

with Q = 10 H 2 O molecules ster s . . The scale 

1/2 

distance Q , ) The electrons are rapidly cooled and 

all but some 10% of the ions undergo charge exchange. 

This behaviour is not sensitive to our assumption (15) 

for the flow velocity, since it occurs explosively 

4 

quickly. We conclude that this 1-2x10 km region is 
effectively a transition region over which the outer 
plasma carrying the energy and ion flux of the solar 
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wind changes continuously to plasma created and 
energised by the solar radiation. The purely 
cometary ionospheric plasma, flowing outwards with the 
expanding gas coma, would have stagnation pressure only 
10% or less of that of the solar wind at the transition 
position — it can hardly affect the flow there. 

Although a stagnation region must occur in the plasma 
flow at some smaller radius , there will be no 
"tangential discontinuity" between plasmas of different 
nature or velocity. 

An important characteristic of the ionospheric 

plasma is that the photo-electrons can cool rapidly to 

\ 

thermal energies before recombining. Rotational 
excitations of H 2 O or OH are effective in the case 
considered. However, if the coma consisted for example 
of pure CO , the cooling mechanism would be more complex 
(section 2), with plasma turbulence trying to thermalize 
an anisotropic distribution of electron energies. The 
corresponding plasma pressure and density might be higher 
and significantly affect the transition flow. But 
in the H 2 O comet, the conclusion is clear, that the 
pressure of the ionospheric plasma is unimportant. 

We have assumed a model coma heated by photo- 
dissociations of H 2 O, this model having a higher 
expansion velocity and temperature and larger ionospheric 
stagnation pressure. If there is no such heating, the 
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plasma pressure would be lower. Shimizu (1974) has 

questioned the reality of the heating in the H^O coma, 

on the grounds that rotational excitations rapidly remove 

the energy of the H-atoms. Indeed, the energy transfer 

from l-2eV H-atoms to the rotational mode appears to be 

comparable to the elastic transfer to translational energy 

is higher by a factor 10, but the energy transferred 

is about 0.025eV rather than 0.2-0. 4eV). This 

indicates that part of the photo-dissociation energy is 

available for heating and increasing the expansion 

velocity of the coma. The conclusion that the coma 

temperature is very low (Shimizu 1974) depends on the 

achievement of thermodynamic equilibrium between the 

rotational levels of H 2 O, and is inapplicable at the 

8 “3 

relevant coma densities (Table 2) of 10 cm , or less. 

The plasma interaction with the coma gas imposes 

strong limits on the place of origin of cometary ions 

which are to form tail rays. For plasma moving at 

4 

around lOkm/s velocxty wxthin the , I =2x10 km scale 

•' X e 

is frictionally decelerated, strongly cooled and liable 
to recombination long before it can flow away. It 

appears impossible for plasma to emerge from inside 

4 . . 

10 km radxus to form taxi rays and streamers. In 
the transition region at 1.5-2. 5km radius, the plasma 
can be cooled to give increased density and still flow 

away before recombination occurs. As such plasma 
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expands adiabatically into tail rays , the recombination 


rate per unit mass changes as 

n n -kCyD+J- 

pa pT p ' , 


(22 ) 


decreasing with p for k ~ j (Table 1). Recombination 
decreases in importance, despite the adiabatic cooling. 

This confirms assumptions of the earlier analysis 
of a tail ray (Wallis 1967) as a jet of plasma, initially 
cold but not undergoing recombination, ejected into the 
solar wind plasma where it is conductively heated and 
frictionally accelerated. The particular transport 
coefficients assumed were based on the transverse 
instabilities of velocity anisotropies in an unmagnetised 
plasma, on which much work has been done recently 
(Davidson 1972). As the magnetic fluctuations were found 
to exceed the expected intensity of any large-scale field, 
the unmagnetised ion-ion instability is indeed 
appropriate, but the postulated electron-ion instability 
may be eliminated by electron gyro-radius effects. The 
order-of-magnitude linear growth rate is, however, 
unchanged. Moreover, the demonstration that the non- 
linear process limiting the instability growth is ion 
trapping (Davidson 1972), confirms the earlier assumption 
(Wallis 1967) equating the growth rate to an effective 
collision or "bounce" frequency. So we consider that the 
earlier results need little modification. They imply, 
we recall, that there was substantial extra mass over and 
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above the observed in the tail rays examined - this 

might well be C"*”, or OH^ . 

+ 4 

Values of the CO density in the coma at 10 km 

- 3 

radius have been given by Arpigny (.1965) as 400cm in 

- 3 

comet Bester Cat l.Oa.u.) and 600-1000cm in comet 

Humason (2.6a.u). These are the same order as the 

ionospheric density (Table 2), although with the lower 

ionization rate at 2.6a.u., the CO production rate would 

29 -1 -1 

have to be higher than 10 ster s . Alternatively the 

transport effects of ions being swept in by the solar 

wind flow can enable higher densities to be reached. It 

is noteworthy that 'envelope' ari^d jet structures were 

4 

observed in comet Bennett at 1-3x10 km ahead of the 

nucleus (Wallis 1973b). The appearance of structures 

at this position corresponds well with the present 

argument that cooling is important in allowing 

condensation of the plasma swept in with the solar wind. 

However, the mechanism for producing structures rather 

than continuous flow is not yet explained. 

When ion-electron recombination is the dominant 

loss process, a recombination instability exists 

(D'Angelo 1967 ) if the coefficient a T varies rapidly, 

- k (y -1) + 2 < - 1 . (23) 

e 

If electron en.ergies were as high as thermal energies 

4 

0.25eV at 10 km (Table 2), the index may be as high as 
k = 2 (Leu et al • 1973), and the plasma might thus be 
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unstable to compressional waves along the magnetic field 


CYg = 3). However, the energy transfer due to rotational 
excitations would exceed that due to recombinations by a 
factor 

6 X 10 ^ e ^ n/s X 10 ^ e^n, 

approximately 500 at 10 km (Table 2). The recombination 
instability might still operate far out in the coma and 
perhaps lead to the formation of 'knots' and other irregularities 
in tail rays. But some other process must underlie the 
formation of envelopes , probably a combination of 
dynamical with ionization and cooling effects. 
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DISCUSSION 

H. Keller ; The outflow velocity used by Wallis is ~ 3-4 km~^ (for H^O, 
OH. . ) on the argument of heating by dissociative excess energie^s whereas 
some observations show only ~lkm s“* (that means no heating of the neutral 
component). Observations are necessary. 

H. U. Schmidt: There may be a misunderstanding. In my discussion I 
assumed a negligible contribution to the temperature from the electrons, so 
that the pressure on the contact surface comes from the expansion of the re- 
maining ions with 1 or 3 km/sec. I think your calciilations are extremely valu- 
able for another purpose, too, i.e. , the electrical conductivity which can be ob- 
tained is important in the same context. 

M. K. Wallis ; Well, I agree that I've ignored things like electron con- 
ductivity. One can take the view that conductivity is high along the field 
lines. I would rather take the view that the plasma is rather turbulent and the 
conductivity on the field lines is not going to be that much different from conduc- 
tivity across the field lines. I agree this is speculation and that it is something 
that needs to be looked into at some stage. 

When you use conservation procedures like this, then you've got to be on 
your guard against that. But the ion pressure, I thought I understood you to say 
earlier that the momentum contribution of the outflowing ions, was unimportant. 
It was more the magnetic stresses which were bigger in effecting the pressure. 

I don't have an outside and inside. There were two models. One is 
flowing in, straight in to the comet and the other one is looking at the plasma 
density in the inner region when you don't have any addition of plasma flow in 
It's just from the photoelectron plasma. 

Now, these two regions have to be matched, of course, and you will have 
some ion pressure. But I'm cooling my electrons down so fast that I'm going 
to recombine the ions. . 

Now, it may be if you add that in it doesn't — that you get a bigger con- 
tribution. I'm not clear on that. We'll have to see. 
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